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ABSTRACT 


A two-dimensional problem of diffraction of a plane 
electromagnetic wave on a smooth 2lt— periodic surface is 

considered. Numerical algorithm, solving this problem is 

developed. , , 

An inverse problem of determination of the shape of 
2'JC-per iodic surface using the performance data of reverse 
scattering is considered. 

Inverse problem was solved by means of minimization of 
the residual functional with the help of gradient descent 
method. The initial data were calculated with the help of 
the numerical method. On each step of iterative method of 
minimization .the residual functional was calculate 
approximately with the help of small slope method. The 
examples of the shape determination are considered. 
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INTRODUCTION 


The new approximate methods, solving the problem of 
diffraction on a smooth two-dimensional infinite wave-1 ike 
surface (for example [1], [2]) give us hope of solving the 
problem of the wavy shape determination using the 
performance data of reverse electromagnetic scattering. 

The aim of this paper is to show the advance of the 
approximate method of small slope in connection with the 
inverse diffraction problem. For this purpose we consider 
only periodic surfaces because for such surfaces there may 
be developed numerical methods, solving the direct problem 
with the high accuracy, so the accuracy of the approximate 
method, solving the inverse problem, may be investigated. 

In the first part of this paper the two-dimensional 
direct problem of diffraction of a plane electromagnetic 
wave on a smooth 2H:-per iodic surface is considered, the 
numerical method, solving the direct problem of diffraction 
is developed: 

With the help of Green's function of Flocke canal the direct 
problem is reduced to the one-dimensional integral 
equation. The kernel of the integral equation contains 
logarithmic singularity, which is expressed in the explicit 
form. The well convergent series for calculating the kernel 
of the integral equation are developed. The integral 
equation is solved using the method of moments . 

In the second part of this paper the small slope method 
is applied for the problem of diffraction on a periodic 
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surface, the algorithm, solving the inverse problem, is 
developed, the examples of shape determination are 
presented: 

The inverse problem of determination of the shape of 
2%-periodic surface using the performance data of reverse 
scattering are considered. Inverse problem was solved by 
means of minimization of the residual functional with the 
help of gradient descent method. The initial data were 
calculated by means of multiple solving of the integral 
equation . On each step of iterative method of minimization 
.the residual functional was calculated approximately with 
the help of small slope method. The formulas for approximate 
calculations of residual functional are presented. 


1, Mathematical formulation of the direct problem. 

The unknown function u satisfies the Helmholtz 

equat i on 


A u 


-.2 

+ ft • u 


0 


( 1 ) 


in the region Q = 1 (x.y) j -oo < 

x. CO 

Here ft is wave number, R=- , 

function . 

The boundary condition for 


x < f (y) . 0 < 

f(y) is smooth 

the function u 


y < 2% > . 
2Tt-per iodic 


is 


u(f (y) ,y) = 0 


L 


In the region 
radiation condition 


x < x Q = 


inf 

[ 0 . 2 %! 


f(y) the 


+oo -ir x 

ift(x- cosa+y ■ sina) T n 

u = e + L x n e 

n=-oo 


iX y 
n 


(3) 


is imposed on u . Here a is the angle between the wave 
vector of incident wave and 

T 


x - axis, X = ft- sina + n. 


r n =h 2 - x^ . Re r n ;> 0 Im r n > 0. 
amplitudes of scattered plane waves. 


n 


are unknown 


The function 


is also assumed to satisfy the 
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Flocke conditions: 


(x, 2%) = u(x. 0) • e 


du du it 

59 (X. 2%) - ^ (x. Ol e 


where t = 2 reft-sina. 

2. Mathematical formulation of the inv erse problem^ 

Let us consider the set of direct problems (1) - (5) 

for resonance values of parameters k and a 


( a e (~.0)u(0.j) . 


fe = = ft(a) = 


Problem: 


For the given function T . =lgnc » “ T ^ lv na ( 1 1 L ar * 

determined in (3) ) where « e (-J.O)u(O.J) determine 


2lt-penodic function f(y). 


ANALYS I S 


1. Numerical algorithm, solving the direct pro blem^ 
With the help of Green's function of Flocke canal 


+00 I 

,P) = ^ H c £ fe J <5x 


2 + ( <5y+21tm ) 2 J e 1 


-imt 


m=-oo 


( here dx - x^Xp. dy - V M -Vp. h“'(x> - is the Hankel 

function of the first kind of order 0 ) the problem is 

reduced to the one-dimensional integral equation for the 

^(f(Y),Y) = 


J 
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du 


^ * y M* ' Y M' ^ y P^ 


• y p]' I( VW fl V-V dY p 


1 


= e 


ift( f (y ) • cosa+y ■ sina] 

M M 


( 8 ) 


-fT 


Here I (y ) = 4 1 + [ f (y ) 

The integral equation (8) was solved with the help of 
me thod-of -moments : 

Let us divide the segment [0,27c] into N equal length 
segments, using points y 
functions : 

1 , y 


( y “0, y =27C ). Consider 

O N 


<P (y) = < 


[ y ,y 1 

i-l X. 


0 


y e [ y ,y ] 

1-1 V 


Let us seek an approximate solution of equation (8) in 
the following form : 

N 


( y ) * J D* * (y 


(9) 


V — 1 


where coefficients D. are to be determine. 

i. 

Function ^ N (y) assume to satisfy equation (8) in points 

1 N 

y * -* ( y + y. ) . It gives following equations for D 

1 2 L-l l. V, 

i — 

2 

coefficients determine : 


N 


y 

j 


2 d " f G [ f(y >>- y .-f'V-yp] i<y p > dy t 

t = 1 J V L L 
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r = exp f ift(f(y ( ) • cosa + y ^sina J (10) 

V l t 

2 2 


The expression (11) gives the well convergent series, which 
gives us the method for calculating the kernel of the 
integral equation (8) . 


iX 6y iy | <5x | 

O O' 1 

G ( M.P ) = — — + 

2x o 

oo i^<5y 

+ ^ [ e b [ ch (~ | <5x | ) * ( n , 6x , 6y ) - 

n=l 

- i • sh(^|«5x| ) -§ 2 (n,6x.6y) J + R(n,6x.6y) j 
Here M = (x || .y |€ ) . P = (x p' y p } ' b=2%. 


§ ( n , 6x , 6y ) 

i 


cos (n<Sy ) -n | <5x | 
n 


(12) 


, , r f , sin(n6y) ^-n|6x| 

$ (n,(5x,i>y) = -e 1 1 

2 n 


(13) 


R(n,6x,6y ) 


-H 


iX n 6y ik n |^x| iX - n 6y lx - n l 6x l 


n 


-n 


)- 


i^6y -n|6x| 

e _e 

n 


ch(^-|<5x| ) - cos (n6y) - 


sh(^-|6x| ) ■ sin(n<5y) J 


R satisfy the expression (15) 


(14) 


L 


J 


568 

Third lntemauona! Conference on Lnverse Design Concepts and Optimization in Engineering Sciences 
(ICIDES-Tin. Editor: G.S. Duiikravich. Washington D.C.. October 23-25. !°91. 


r 




| R ( n . 6x ,6y) | < - 


n 


( 15 ) 


Rows (12), (13) can be summed: 


00 


^ ( n , 6x , 6y ) = - 


1 n2 + 


I <5x| 


n-1 


2 . 6 y 
+ sin ( 


% > ] 




,2 / , 

sh ( — ) + 


( 16) 


00 


J $ z (n,6x,6y) = arctg 


sin (6y ) 


n=l 


<5x \ 

e 1 1 - cos (6y ) 


(17) 


The kernel of the integral equation contains 

logarithmic singularity, which is expressed in the explicit 
form in ( 16 ) . 

2. Approximate method, solving direct problem 

The well known small perturbation method gives 

following formulas for amplitudes T : 

n 

T S ‘ P ‘ « A 6 + B f ■ + F C -f f (18) 

n on n n ** mn m n-m 

m 


Here f - — 
n 21U 


2 % 

J e inY f(y) dy 


. * -f l 
on 1 0 


n = 0 
0 n * 0 


A 1, B = - 2i r 

n o 


C =2 r r 

mn o m. 


The small slope method solving the problem of diffraction 
on a smooth two-dimensional infinite wave-like surface was 
presented in [1]. In the case of periodic surface amplitudes 


n 


were sought in the following form: 


f S.S. 
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r 


i 


(a + £ b -f-e 1Y )dy 
n ^ mn m 
m 


(19) 


where a . b were constants to be determined 
n mn 

s s 

values T ' ' satisfy the following conditions: 

n 

a) |t S ' S ‘ - t S ' P '| < const-*: 3 , c - 0 
1 n n 1 

where e = max | f (y) | 

0<y<21t 


Let the 


( 20 ) 


b) Shift surface along x - axis condition: 
S.S. 


T“ ‘ ‘ [ f ( y ) +x ] = e 

n o 


1 ^ r n +r r )X o q s 

O r, °_ T S.S. [f (y) ] 
n 


( 21 ) 


c) shift surface along y - axis condition: 

C q in y« q q 

T S ' S ' [ f (y+y )] = e -T 3 ' b [f(y)] (22) 

no n 

Condition (20) will be valid for (18) if 

a = A (23) 

o 


i(r + y ) • a + b = B 
o x\ n nn n 


(24) 


^mn ^n-m , n 


(25) 


where Q mn - C mn 


^ — < 2l b mn " ( ^0 + >V a n >* and the 


conditions (12) and (13) will in turn be satisfied if 


b = 0 
on 


(26) 


In the article [1] the following expressions were proposed. 
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which satisfy (23) - (26): 


a 


n 




c (C + C 

2 on nn 


2 1 (y + 
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o n 
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mn 
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2 tr n +r 0 ] 
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n-m , n 
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+ C 
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3. Algorithm, solving the inverse problem . 
Let 


d = (d , • • * , d 

1 M 


, d ) 

2M + 1 


Consider the set of surfaces , which are determined with 
the help of functions: 

m m 

f (y ) = ) d -sin ny + d + ) d -cos ny 

L n M + l L n+M+1 

n= 1 rt~ i 

Let us suppose that we know values of the inverse scattering 
function for surface, corresponding to vector d° : 


T° = T ( (X ) , T° = T«X ) , 

1 12 2 


• • • , T° - T(d ) 

N N 


for angles d i ,d 2 , 


d 


N 


The numerical algorithm for vector 


reconstruction 


on values T° , ••• T° was constructed. This algorithm is 

1 N 

based on minimization of the residual functional with the 
help of gradient descent method. Direct problem on each 
iteration step is solved with the help of small slope 
approximate method. On each step of iterative method of 
minimization , the residual functional was calculated 
approximately with the help of small slope method. 

Residual functional is determined by following 
expression : 
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-y -> 

$ (a.d) 


N 

1 

m= i 


t - Rta 

m rr 


.d) | 


(27) 


1 


Here & = ( . • Ct N ) , 

R(a.d) is the inverse scattering functional 

corresponding to small slope method: 


R(a 


,d) I 2 % ( 6 OTn 


I ctg a| 

+ if • ( - r ) 

m rn 


) 


m 


2 7T 


J 


i(m+signa)y i|ctgCL|f(y) 
e • e dy 


277 

*- - zft-J . 

o 

Coefficients T° were calculated with the help of numerical 

m 

method ( 7 ) - ( 17) . 

RESULTS 

Consider the example of solving the direct problem with 
the help of IBM PC AT 386-387 computer. The values of 

parameters are : ft = 1.4 , Ct = 30° , f(y) = 0.3- sin y ; 
in numerical algorithm N = 50 . 


Numerical Small slope 

algorithm method 


energy error 

1 

o 

* 

CO 

o 


0.16-10 2 

CPU time 

3 min 54 

sec 

1 sec 

T 

-2 

-0.072 


-0.073 

T 

-i 

0.349 - 

0 . OlOi 

0.347 - 0 . 005 i 

T 

o 

-0.926 + 

0 . 052i 

-0.926 + 0 . 049 i 
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^ Consider examples of Furie coefficients reconstruction 

of surface determine function: 

-» 

First simulation was carried out for N=10 , M=1 ,d 

■** 

( 0, -0.3, 0.51 ) . After 7 iterations vector d : ( 0.000, 
-0.300 , 0.503 ) was obtained as the minimum of the 

functional (27) . 

Second simulation was carried out for N=6 , M=2 

-» 

d° = (0, 0, -0.1, 0.2 0.2 ) . After 44 iterations vector 

-** 

d : ( 0.000, 0.000, -0.092, 0.195. 0.190 ) was obtained 

as the minimum of the functional (27). 


REFERENCES 

1. A . G . VORONOVICH, Small slope approximation in the 

theory of scattering of waves at rough surfaces . 
Journal experim. and theor. Phys . , V. 89(1985), issue 

1(7), pp . 116— 125(In Russian.) 

2. Ishimary A., Winebrenner D. Investigation of a surface 

field phase perturbation technique for scattering from rough 
surf aces . //Radio science -1985.- V. 20. -N . 2. pp . 161 - 

170 . 


L 


J 


